Nucleon Mass Splitting at Finite Isospin Chemical Potential 
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We investigate nucleon mass splitting at finite isospin chemical potential in the frame of two 
flavor Nambu-Jona-Lasinio model. It is analytically proved that, in the phase with explicit isospin 
symmetry breaking the proton mass decreases and the neutron mass increases linearly in the isospin 
chemical potential. 
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Recently, the study on QCD phase structure is ex- 
tended to finite isospin density. The physical motivation 
to study QCD at finite isospin density and the corre- 
sponding pion superfluidity is related to the investigation 
of compact stars, isospin asymmetric nuclear matter and 
heavy ion collisions at intermediate energies. While there 
is not yet precise lattice result at finite baryon density, it 
is in principle no problem to do lattice simulation at fi- 
nite isospin density [![. It is found thatQ the phase tran- 
sition from normal phase to pion superfluidity happens 
at a critical isospin chemical potential which is exactly 
the pion mass in the vacuum, /ij = m^. The QCD phase 
structure at finite isospin density is also investigated in 
many low energy_efFective_models, such as chiral pertur- 
bation theor 
model 



HHHl, Nambu-Jona-Lasinio (NJL) 
random matrix methodjllj, ladder 
QCD 12], and strong coupling lattice QCD|13j. 

There are two kinds of isospin symmetry breaking. At 
low isospin chemical potential, /i/ < /ij, the isospin sym- 
metry is explicitly broken. In this phase, while a and 7r 
mesons are still the eigen collective modes of the Hamil- 
tonian of the system, the n meson splits into three chan- 
nels with linear fij dependence for n + and n_ masses [fOj. 
M 7r+ = m n — fii and M w _ = + fj,j. At high isospin 
chemical potential, /i/ > fj, c Iy the isospin symmetry is 
spontaneously broken with nonzero pion condensate. In 
this phase, a and 7r are no longer the eigen collective 
modes of the system, and the new eigen modes are the 
linear combination of them. In the pion supcrfluid, the 
simple linear dependence of the collective modes on the 
isospin chemical potential disappears, and masses of the 
new eigen modes are determined by the complicated pole 
equations [lOj. 

Like the three 7r mesons which form an isospin triplet, 
proton and neutron are different isospin states of nu- 
cleon. Similar to the linear fii dependence for charged 
pions, we guess that there may exist a simple /i/ depen- 
dence for proton and neutron masses in the phase with 
only explicit isospin symmetry breaking. In this paper, 
we calculate the nucleon mass splitting at finite isospin 
chemical potential in the NJL model. We first extend 
the nucleon structure [3] of a quark and a diquark in 
the model to dense medium, and then determine analyti- 
cally the proton and neutron masses at finite baryon and 
isospin chemical potentials. 



We take the two flavor NJL model defined by the La- 
grangian density 
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where C = 17270 is the charge conjugation matrix, 
the superscript t denotes the transposition operator, 
niQ is the current quark mass, /x is the quark chemi- 
cal potential matrix in flavor space /i = diag (fi u , ji^) = 
diag (/x.b/3 + /Ar/2, /is/3 — /i//2) with and fii being 
baryon chemical potential and isospin chemical potential, 
G s and Gd are coupling constants in color singlet chan- 
nel and anti-triplet channel, r = (ti,T2,T3) are Pauli 
matrices in flavor space, and ej and e c are totally anti- 
symmetric tensors in flavor and color spaces. 

Since we are interested in the nucleon mass splitting in- 
duced by explicit isospin symmetry breaking at fij < /ij, 
there are no pion superfluid and color superconductor, 
the only condensate is a = (tptp) which is the order pa- 
rameter of chiral phase transition. 

In mean field approximation, the thermodynamics of 
the system can be evaluated asflfjf 
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in terms of the quasi-particle energies Ef = E p ±/ii with 
E p = -\/p 2 + to 2 , where m = too — 2G s a is the effective 
quark mass. From the thermodynamic stable condition 
of the system, the physical condensate is determined by 
the gap equation 



on 







(3) 



which gives the temperature and chemical potential de- 
pendence of the quark mass, to(T, fij). At zero 
temperature, the quark mass keeps its vacuum value 
in the normal pha se without pion superfluid and color 
superconductor [lOj, to(0, fii) — m(0,0,0). 

Since we still have color symmetry, the effective quark 
propagator S in mean field approximation is independent 
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of color index and diagonal in flavor space, 

S (p, p s , M/) = dia 9 (Su Mu) i (p, Md)) (4) 
with the elements 

= S^ 1 (po + Mi, P>0) , i = u,d. (5) 

In the NJL model, the meson modes and diquark 
modes are regarded as quantum fluctuations above the 
mean field, they can be calculated in the frame of 
rpa0, EE E3. When the mean field quark propaga- 
tor is diagonal in flavor and color spaces, e.g. the case 
with only chiral condensate, the summation of bubbles 
in RPA selects its specific channel by choosing at each 
stage the same proper polarization function, the diquark 
propagator can be written as 



Dd (k,(J,B,(J>l) = 



1 - 2G d ILj(fc,MB,M/) 
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where LT^ is the diquark polarization function constructed 
by two quarks [14j. 
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There are three kinds of diquarks dj, <% an d <% con ~ 
structed by quarks with colors 2 and 3, 1 and 3, and 
1 and 2, respectively. Since the quark propagator is color 
independent, they are degenerated. Taking the trace in 
flavor and color spaces and employing the relation 



CS(p,fiB,fJ>iYC 1 = S(-p,-/j, B ,-Hi), 
the polarization function can be simplified as 

d 4 p 
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now the trace is taken only in Dirac space. 

Taking the replacements po ~ > Po — M« f° r t ne first term 
in the square bracket and po — > po — Md f° r the second 
term, the chemical potential dependence of lid and D d 
can be reflected in a shift of the diquark energy, 



Ud(k,fJ,B,fJii) = n d (fc + -MB,k, 0, 0), 
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D d (k,HB,m) = D d (k + -Ms, k, 0,0). (10) 



It is easy to see that the diquark polarization and prop- 
agator are isospin chemical potential independent. This 



is easy to be understood. From the interaction in di- 
quark channel in (fT]), only the quarks with different fla- 
vor and color can construct a diquark, and a u quark 
and a d quark have opposite isospin chemical potential. 
We should emphasize two points for the relation (| 10[) . 
First, (fTU|) is valid only in normal phase where the effec- 
tive quark mass m is a constant. In the pion superfluid 
phase, m is a function of chemical potentials, m(ps, pj), 
and then the relation (fTTjf fails. Second, at finite tem- 
perature the Lorentz invariance is broken down, and the 
integration over quark energy po or diquark energy fco 
is changed into a summation over the discrete Fermi or 
Bose frequencies. Therefore, the replacement used to get 
dTTJj) does not work at finite temperature. 

Now we investigate nucleons in the NJL model. A 
nucleon constructed by three quarks can be described 
by the three-body Faddeev equation [li]]. In the static 
approximation [1 7j . the Faddeev equation is reduced to 
an effective BS equation constructed by a quark and a di- 
quark, the quark propagator between two quark-diquark 
bubbles becomes a constant — 1/m, and then we can use 
RPA again to derive the nucleon propagator [l8j 
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where II is the nucleon polarization function 
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xSi(q -k + Mi, q-k,0). (12) 

Since the diquark contains a u and a d quark, we take 
Si = S u for protons and Si = S d for neutrons. Making a 
transformation for the diquark energy, ko — > fco — 2p_e/3, 
we have 



IIjv(g, Mb, M/) 
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xSi(q - fc + -Mb + p/, q - k, 0). 



The chemical potential dependence of LT^ and is only 
reflected in the quark propagator. 

From the definition of particle mass, nucleon mass Mm 
corresponds to the pole of the nucleon propagator at zero 
momentum, 



1 H U N (q,fj, B ,Hi) 

m q =M N ,q=0 

In the vacuum it is reduced to 



0. 
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for the vacuum mass mpf. At this stage, we can deter- 
mine the four parameters in the NJL model, the current 
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quark mass mo, the momentum cutoff A and the two cou- 
pling constants. By fitting the pion mass m v = 140 MeV, 
pion decay constant f v = 93 MeV and the quark mass 
to = 400 MeV in the vacuum, we have too = 5 MeV, 
A = 593 MeV and G s = 6.92 GeV~ 2 . The coupling 
constant Gd in the diquark channel is controlled by the 
nucleon mass in the vacuum. The ratio between the two 
couplings is Gd/G s — 1.1, corresponding to to at — 940 
MeV and diquark mass Md = 600 MeV which can de ob- 
tained from the pole of the diquark propagator © . Note 
that the value of the ratio is in between the estimated 
minimum and maximum values [l6| . 

From the comparison between the pole equations (fl4|) 
in matter and (| 1 5(1 in vacuum, we obtain analytically the 
relation between the nucleon masses Mn in matter and 
tojv in vacuum, 

2 

+ -fiB + m = Wni i = u,d. (16) 
o 

For a proton which is a u(ud) state and a neutron which 
is a d(ud) state, we have 

M p = TO TV - /is - -/!/, 

M n = to 7v - Ms + ■ (17) 

Just as we expected, the nucleon mass splits into two 
channels in matter at finite isospin chemical potentia 
/ij < /Xj, and the proton and neutron masses change 



linearly in /i/. Again, we should note that these linear 
relations are true only at zero temperature and in normal 
phase with only chiral condensate. 

In summary, we have investigated the nucleon mass 
splitting in the phase with explicit isospin symmetry 
breaking in the frame of two flavor NJL model at zero 
temperature. While proton and neutron arc still degener- 
ated at finite baryon chemical potential hb and the mass 
drops down linearly in fig, the introduction of isospin 
chemical potential fii leads to their mass splitting. The 
proton mass decreases and neutron mass increases lin- 
early in fii. This remarkable mass splitting may sig- 
nificantly change the properties, such as the threshold 
values, of those dynamic processes including protons and 
neutrons in isospin asymmetric nuclear matter. 

At finite temperature, the pole equation for proton and 
neutron masses is still valid but it needs numerical calcu- 
lations to obtain their temperature dependence. When 
we go into the pion superfluid region with finite pion 
condensate, the nucleon mass splitting is still true, but 
the linear relation fails, because the quark propagator is 
off-diagonal in flavor space and the diquark propagator 
and in turn the nucleon propagator become much more 
complicated. 
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